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Abstract. In recent work on nonequilibrium statistical physics, a certain Markovian exclusion model called an asym- 
metric annihilation process was studied by Ayyer and Mallick. In it they gave a precise conjecture for the eigenvalues 
(along with the multiplicities) of the transition matrix. They further conjectured that to each eigenvalue, there cor- 
responds only one eigenvector. We prove the first of these conjectures by generalizing the original Markov matrix 
by introducing extra parameters, explicitly calculating its eigenvalues, and showing that the new matrix reduces to 
the original one by a suitable specialization. In addition, we outline a derivation of the partition function in the 
generalized model, which also reduces to the one obtained by Ayyer and Mallick in the original model. 

Resume. Dans un travail recent sur la physique statistique hors equilibre, un certain modele d' exclusion Markovien 
appele "processus d' annihilation asymetriques" a ete etudie par Ayyer et Mallick. Dans ce document, ils ont donne 
une conjecture precise pour les valeurs propres (avec les multiplicites) de la matrice stochastique. Ils ont en outre 
suppose que, pour chaque valeur propre, correspond un seul vecteur propre. Nous prouvons la premiere de ces 
conjectures en generalisant la matrice originale de Markov par l'introduction de parametres supplementaires, calculant 
explicitement ses valeurs propres, et en montrant que la nouvelle matrice se reduit a l'originale par une specialisation 
appropriee. En outre, nous presentons un calcul de la fonction de partition dans le modele generalise, ce qui reduit 
egalement a celle obtenue par Ayyer et Mallick dans le modele original. 

Keywords: Reaction diffusion process, non-equilibrium lattice model, transfer matrix Ansatz, partition function, 
characteristic polynomial, Hadamard transform. 



1 Introduction 

In the past few years, special stochastic models motivated by nonequilibrium statistical mechanics have 
motivated several combinatorial problems. The most widely studied problem among these has been the 
totally asymmetric simple exclusion process (TASEP). The model is defined on a one dimensional lattice 
of L sites, each site of which either contains a particle or not. Particles in the interior try to jump with 
rate 1 to the site to the right. The jump succeeds if that site is empty and fails if not. On the boundary, 
particles enter with rate a on the first site if it is empty and leave from the last site with rate (3. This was 
first solved in 1993 by developing a new technique now called the matrix product representation [1 1. 

It was initially studied in a combinatorial setting by Shapiro and Zeilberger in an almost forgotten 
paper [3| in 1982, but only after the steady state distribution of the model was explicitly presented in 
IH, the problem gained widespread attention. One of the reasons for this interest was that the common 
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denominator of the steady state probabilities for a system of size L was Cl+i, the (L + l)-th Catalan 
number. One of the first articles to explain this fact combinatorially was the one by Duchi and Schaeffer 
101, who enlarged the space of configurations to one in bijection with bicolored Motzkin paths and showed 
that the steady state distribution was uniform on this space. The analogous construction for the partially 
asymmetric version of the model (PASEP) has been done in |5 1. 

Further work has been on the relationship of the total and partially asymmetric exclusion processes to 
different kinds of tableaux by Corteel and Williams (6J [7j [8) (permutation tableaux, staircase tableaux) 
and by Viennot [9| (Catalan tableaux), to lattice paths [10], and to Askey-Wilson polynomials [7|. 

Just like the common denominator for the TASEP of size L was the Catalan number Cl (which has 
many combinatorial interpretations), the common denominator for the asymmetric annihilation process 
considered in [2 1 in a system of size L at a = 1/2, /3 = 1 is 2^ 2 / which is the number of domino tilings 
of an Aztec diamond of size L as well as the number of 2-enumerated L x L alternating sign matrices. 
One can therefore hope to enlarge the configuration space as was done for the TASEP |4] to explain this 
phenomena. 

The remainder of this extended abstract is organized as follows: In Sec. |2]we describe the model of the 
asymmetric annihilation process. In Sec. [3]we present some of the main results obtained by Ayyer and 
Mallick in [2]. Their work lead to a conjecture about the spectrum of this process. In Sec. |4]we prove this 
conjecture by appropriately extending the model and viewing it in a different basis obtained by a variant 
of the Hadamard transform. In the concluding section we outline the derivation of the partition function 
for the generalized model using the same transformation, an approach very different from the way Ayyer 
and Mallick obtained the partition function in the original model. 



2 The model 

Motivated by Glauber dynamics of the Ising model, Ayyer and Mallick [2| considered a non-equilibrium 
system on a finite lattice with L sites labelled from 1 to L. States of the system are encoded by bitvectors 
b = b\bi ■ ■ ■ of length L, where bj € B = {0, 1}, so that we have a total of 2 L states. These bit vectors 
may be represented numerically using the binary expansion (6)2 = 2 1 +&l_2 2 2 + - • -+b\ 2 i ~ 1 , 

which introduces a total order on B L , so that we shall write b < c iff (6)2 < (0)2. All matrices and vectors 
are indexed w.r.t. this order. 

The evolution rules of the system introduced in [2 1 can now be stated as rewrite rules for bit vectors: 



• In the bulk we have right shift and annihilation given by 

right shift 10 01 with rate 1, 
annihilation 11 — > 00 with rate A, 

and visualized (for L = 8) in Fig. [T] 



Fig. 1: Right shift 00110101 -» 00110011 and annihilation 00110101 -» 00000101 
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• On the left boundary, particles enter by left creation in a way consistent with the bulk dynamics. A 
particle at site 1 may also be left annihilated (due to a virtual particle at site 0). Therefore, the first 
site evolves as 

left creation — >• 1 with rate a, 
left annihilation 1^0 with rate aX, 

as illustrated Fig. [2] 

* >• -O »» » • -O »• > 0— >•— <! O— -O— O ••— O ■> 'O ■> 

Fig. 2: Left creation 00110101 -> 10110011 and left annihilation 10110101 00110101 

• Particles can exit from the last site by right annihilation (with a virtual particle at site L + 1) 
according to 

right annihilation 1 — > with rate /?, 

as illustrated by Fig. [3] 

Fig. 3: Right annihilation 00110101 -> 10110010 



Note that all transition rules except left creation are monotonically decreasing w.r.t. the natural order of 
bit vectors. Thus the transition matrix, as discussed in the next section, is not in triangular shape. 
Following [2], we will take A = 1 as that is the only case for which they derive explicit formulae. 

3 Algebraic properties of the model 

Is this section we present without proofs the main results as obtained by Ayyer and Mallick in [2|. First 
recall the general concept: 

Definition 1 A (continuous-time) transition matrix or Markov matrix or stochastic matrix is a square 
matrix of size equal to the cardinality of the configuration space whose (i,j)-th entry is given by the rate 
of the transition from configuration j to configuration i, when i is not equal to j. The (i, i)-th entry is then 
fixed by demanding that the entries in each column sum to zero. 

The Markov chain we defined in the previous section satisfies what [2 1 call the "transfer matrix Ansatz". 
The following general definition applies to any family of Markov processes defined by Markov matrices 
{Ml} of increasing sizes (in most physical applications, L is the size of the system). 

Definition 2 A family Ml of Markov processes satisfies the Transfer Matrix Ansatz if there exist matrices 
Th.L+ifor all sizes L such that 

(TMA) M L+1 T L>L+1 = T L}L+1 M L . 
We also impose that this equality is nontrivial in the sense that Ml,+\Tl l+1 7^ 0. 
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The rectangular transfer matrices Tl,l+i can be interpreted as semi-similarity transformations connecting 
Markov matrices of different sizes. 

The last condition is important because there is always a trivial solution whenever we are guaranteed a 
unique Perron-Frobenius eigenvector for all transition matrices Ml. If \vl) is this eigenvector of Ml and 
(1l\ = (1, 1, ■ ■ ■ , 1), the matrix Vl,l+i = \ v l+i)(^-l \ satisfies (TMA) since the Markov matrices satisfy 
the conditions (1l\Ml = and M l+\\v l+i) = 0. 

The above definition leads immediately to a recursive computation of the steady state vector, which is 
the zero eigenvector. First we have 

= Tl,l+iM l \vl) = Ml + iT l ,l+i\vl), 

which, assuming Tl.l+i\ v l) 0, and taking into account the uniqueness of the steady state, allows us 
to define so that 

T l ,l+i\vl) = \vl+i)- 

This is very analogous to the matrix product representation of [ 1 1 because the steady state probability of 
any configuration of length L + 1 is expressed as a linear combination of those of length L. The transfer 
matrix Ansatz is a stronger requirement than the matrix product representation in the sense that not every 
system which admits the representation satisfies the Ansatz. For example, the only solution for (TMA) in 
the case of the TASEP is the trivial one. 

For our system introduced above, the Markov matrices Ml are of size 2 L . As mentioned, the entries of 
these matrices are indexed w.r.t. the naturally ordered basis of binary vectors of length L. For convenience, 
here are the first three of these matrices: 



Mi = 



-a 
a 



a + /3 
-a-p 



M 2 = 





P 


a 


1 







1 


a 


a 







P 





a 





* 



, M 3 = 



★ 


p 





1 


a 





1 








★ 


1 








a 





1 








* 


p 


1 





a 



















1 





a 


a 











★ 


P 





1 





Q 











■k 


1 











Q 











* 


P 











o 














As for the diagonal elements *, they have to be set such that the column sums vanish. 

We now state without proof some important results on the Markov matrices of the system. These are 
proved in 0. We first show that the Markov matrix itself satisfies a recursion of order one. 



Theorem 1 Let a denote the matrix a 



1 

1 



, and 1 l denote the identity matrix of size 2 L . Then 



Mr 



M L -i - a(a ® l L - 2 ) 


al L -i + (<t®1l-2) 


al L -i 


M L -i - 1l-i - a{<j® 1 L _ 2 ) _ 



where Ml is written as a 2 x 2 block matrix with each block made up of matrices of size 2 



L-l 



The transfer matrices can also be explicitly constructed by a recursion of order one. 
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Theorem 2 There exist transfer matrices for the model. If one writes the transfer matrix from size 2 



L-l 



to size 2 by a block decomposition of matrices of size 2 



L-l 



x 2 



L-l 



[ L-1,L 



(L-l) 



T, 



(£) 



1? 



r. 



(L-l) 



f/ien f/ie matrix Tl,L+i can be written as Tl.l+i 



(L) 



T 



, with 



-lrp(L-l) 



a~ T 



(L-l) 



2T 



(L-l) 



-WL-1) 



a~ L T 



OL~ T, 



T, 



(i) 



2T _(L-1) 


T 2 (L V® l£_ 2 ) ' 





T (L-i) 



This, along with the initial condition 



Ti, 



1 + (3 + aft a + (3 + a(3 
a 1 
a + a/3 a/3 
a 



determines recursively a family of transfer matrices for the matrices M^. 

We can also use the transfer matrices to calculate properties of the steady state distribution of the 
Markov process. One quantity of interest is the so called normalization factor or partition function. 

Definition 3 Let the entries of the kernel \ vl) of Ml be normalized so that their sum is 1 and each entry 
written in rationally reduced form. Then the partition function Z^for the system of size L is the least 
common multiple of the denominators of the entries of \ vl )• 

Because of the way the transfer matrix has been constructed Zl is the sum of the entries in vl. For 

~a + 



example, the system of size one has | v\ ) = 



a 



, whence Z\ = 2a + (3. 



Corollary 3 The partition function of the system of size L is given by 

Z L = 2(V)(1 + 2a) L - 1 (l + p) L - 1 {2a + (3). 

4 Spectrum of the Markov matrices 

In this section, we consider the eigenvalues of the Markov matrices Ml of the asymmetric annihilation 
process. The following result was stated as a conjecture in O . This will be a corollary of the main result 
(Theorem[8]l of this article. 



Theorem 4 Let the polynomials Al(x) and Bl(x) be defined as 
A L (x) 



[L/2] [L/2] 

Jl (x + 2k)( L ™ 1 ) 1 B L {x)= Jl (x + 2k+ l^Vi). 

fe=0 fc=0 



Then the characteristic polynomial Pl{x) of Ml is given by 

P L {x) = A L (x)A L (x + 2a + (3)B L {x + f3)B L (x + 2a), 
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and successive ratios of characteristic polynomials are given by 
Pl+i(x) 



Pl{x) 



B L {x + l)B L (x + 2a + p + l)A L (x + /3 + l)A L (x + 2a + 1). 



This gives only 2L distinct eigenvalues out of a possible 2 L . There is therefore the question of diago- 
nalizability of the Markov matrix. Ayyer and Mallick (2) further conjecture the following. 

Conjecture 1 The matrix Ml is maximally degenerate in the sense that it has exactly 2L eigenvectors. 

For L > 1 we regard B L as the vector space of bitvectors of length L (over the binary field) . The usual 
scalar product of vectors b, c € B L will be denoted by b ■ c. We will take the set Vl = { \ b ) ; b€ M L } as 
the standard basis of a 2 L -dimensional (real or complex) vector space, which we denote by Vl- Indeed, 
we will consider Vl as a vector space over an extension over the real or complex field which contains all 
the variables that we introduce below. To be precise, we take Vl as a vector space over a field of rational 
functions which extends the real or complex field. 

The following definitions of linear transformations, when considered as matrices, refer to this basis, if 
not stated otherwise. Vl is the L-th tensor power of the 2-dimensional space Vi in an obvious way. 



The transformation a of V\ is given by the matrix a = 
mations a b of V L for b = bib 2 ■ ■ ■ b L € B L : 



1 

1 



and this extends naturally to transfer- 



er 6 = f7 blb ^- bL = a bl (g>a b2 ® • •• ® a bh . 
Definition 4 For a vector a. — {a.b)bea L of variables we define the transformation Al{ol) ofVL as 

A L (a) = a bO h . 

b£M L 

A direct way to define these matrices is (b\AL |c) = ab® c , (b, c G B L ), where © denotes the 
component wise mod-2-addition (exor) of bit vectors. 
For 1 < j < L we define the involutive mappings 

</>j ■ : M L -> B L : b x ■ ■ ■ bj-ibj b j+1 . . . b L i-> (j>jb = bi . . . bj-ib]b j+1 ...b L 

by complementing the j-th component, and involutions 

ipj : B L -> B L : b i-> (pjfa+ib 

by complementing components indexed j and j + 1, where ipL is the same as 4>l- 
Definition 5 1. For 1 < j < L we define the projection operators Vl.j acting on Vl by 

VL,i= E \b)(b\-\//(b))(b\ 

beB L 

2. For a vector b — [fii, fa, ■ ■ ■ , Pl) of variables we put Bl{0) = Z~2i<j<L Pi ^L.j- 
E>l (/3) denotes the matrix representing Bl (/3) in the standard basis Vl- 
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Note that in the sum for V L j only summands for which bj = 1, i.e., for which tpj(b) < b, occur. 
Indeed: this condition allows only for two situations to contribute: 



bj bj + i = 10 i y bj bj + i = 01 
bj bj + i = 11 i — ^ bj bj + i = 00 



(right shift) 
(annihilation) 



(1) 
(2) 



Thus these operators i[>j encode the transitions of our model. Also note that by its very definition Bl{(3) 
is an upper triangular matrix. 

Writing (3 = (/3, 7, S) instead of (fa, fa, fa) we have for L = 3 



B 3 (fal,6) 



"0 


-5 





-7 








-0 








5 


-7 














-0 








7 


-5 


-0 




















(5 + 7 





-0 























-5 





-7 

















5 + p 


-7 























7 + /3 


-d 























5 + 7 + /3 



Our main concern is now with the transformation given by 

M L (a,(3)=A L (a)-B L (f3). 

Before we can state the main result we have to introduce some more notation, But before doing so, we 
note that the corresponding matrix M L (a, (3) = A L (a) — B L ((3) reduces to the matrix M L above when 
properly specialized: 

Lemma 5 We have A L (a') — B L ((3') = M L for a' = (a' b ) beM L and (3' = {P'j)i<j<L given by 



a h 



-a if b = 00... 00 

a if b= 10... 00 and ft 

otherwise 



1 ifl<j<L 
ifj = L 



We will now consider the transformation A4l{ol, (3) in a different basis of Vl- Let H = 



be the familiar Hadamard matrix and define Hl as its L-th tensor power, the matrixL Hadamard transform 
of order L: 

1 



H L =H 



2 L / 2 



[(-1) 



be" 



l b.cei 



The columns of this matrix, denoted by | w b ) = \ b ) for b e B L , form an orthonormal basis Wl = 
{Hl \ b) ; b £ B L } of Vl. The following assertion is easily checked: 

Lemma 6 The (pairwise commuting) transformations a c (c £ M L ) diagonalize in the W^-basis. More 
precisely: 

a c \w b ) = {-l) b - c \w b ) (b,c£M L ). 
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Thus also the transformation Al diagonalizes in the W^-basis and its eigenvalues are given by 

(H L -A L -H L )\w b )=\ b \w») 

where \ b = £ ceB i, a c (-l) bc = £ ceB z. a c (b\H\c). 

The crucial observation is now the following: even though the transformation Al diagonalizes in the 
Wl -basis, the transformation Bl doesn't, it is not even triangular in this basis. But it turns out that a slight 
modification of the -basis will be suitable for at the same time diagonalizing Al and bringing the Bl 
in (lower) triangular form. For that purpose we introduce the invertible linear transformation 



A 



. br ^ = 



-> B L : b = 61&2 
where the sum has to be taken in the binary field. As an example (L 



l<i<L-j+l 



l<j<L 



3): 



b 


000 


001 


010 


011 


100 


101 


110 


111 


b A 


000 


100 


110 


010 


111 


011 


001 


101 



The basis Wl = { \ w b )} is nothing but a rearrangement of the Wl -basis, hence the transformation 
Al diagonalizes in this basis as well (with the corresponding eigenvalues). We will write Hl for the 
rearrangement of the Hadamard matrix in this new ordering of the elements of B L . The clue is now 
contained in the following proposition: 

Proposition 7 

H L -B L {f3)-H L = Bl{l3™) 
where (3 rev = (/?l, Pl-X, ■ ■ ■ , Px) is the reverse of (3 — (Pi, P%, . . . , Pl), t denoting transposition. 

We illustrate this proposition in the case L = 3 by displaying matrices 2 Z / 2 H^ (left) and H3 B3 (/?, 7, 8) ■ 
H$ = B^,(8 1 7, (right). Note that Hl is symmetric because A (as a matrix) is symmetric. 



1 1 


1 


1 


1 


1 


1 


1 " 

























1 1 


1 


1 


-1 


-1 


-1 


-1 




-P 


p 

















1 1 


-1 


-1 


-1 


-1 


1 


1 







-7 


7 














1 1 


-1 


-1 


1 


1 


-1 


-1 




-7 





-P 


p + 1 











1 -1 


-1 


1 


-1 


1 


1 


-1 










-8 





8 








1 -1 


-1 


1 


1 


-1 


-1 


1 













-8 


-P 


p + s 





1 -1 


1 


-1 


1 


-1 


1 


-1 




-S 














-7 


7 + 5 


1 -1 


1 


-1 


-1 


1 


-1 


1 _ 







-8 








-7 





-P /3 + 7 + 5 


The proof of Proposition 


[7] will be 


given below. It leads to the 


main 


result by lookin 


g at the matrix 



representation of A^L(a,/3) = Al — Bl in the Wl -basis, where it takes lower triangular form. Hence 
the eigenvalues, which are 2 L pairwise distinct linear polynomials in the a- and /3-variables, can be read 
directly from the main diagonal. In contrast to Theorem|4] the ex-conjecture, all eigenvalues are simple. 
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detM L (a,p) = det[A L (a)-B L {p)]= ]J (V - ' b) 



Illustration of the Theorem for L = 3 recalling f3 = (/3, 7, S) = (/3i, /?2, Ps)'- 



b 


6 A 


AfcA 




(S,%0)-b 


000 


000 1 


+ + + + + + ++] 


■ a. 







001 


100 


+ + + H ] 


■ a 


P 




010 


110 


+ H h+j 


■ a. 


7 




011 


010 


+ H M ] 


■ a 


p- 


H7 


100 


111 


H H h +-] 


■ a 


5 




101 


011 


H + H h] 


■ a 


p- 


V5 


110 


001 


H 1 1 I— ] 


■ a 


7H 


-5 


111 


101 j 


H 1 1 h] 


■ a 


P- 


h7 + <5 



So, as an example, the line for b = 101 contributes the factor aooo — c^ooi ~ a oio — a on + ^ioo — 
aioi — Qfno — am — j3 — 5 to the product. 

To prepare for the proof of Proposition[7]we state without proof simple relations between the transfor- 
mations ipj, <pL-j + i and A: 

Lemma 9 For b,c£M L and 1 < j < L we have 

1. (^-6) A = <j> L -j + i{b A ) 

2. b A ■ ipjc — b A ■ c + &L- 3 +i 
Fact 2. is a consequence of fact 1. 

Proof of Proposition [7] The actions of the transformations Vlj, seen in the Wl -basis, are given by: 



V Ltj :\w b )^ 



' -\ w (4>L- ]+ ib) A } if 6x,-j+i = 0, 
\w bA ) i£b L - j+1 = l. 



To see this, we compute 

-P L ,j\w bA )= ^{b A \H\c)V L , j \c) = (b A \H\c)(\c)-\^c)) 

c£« L c>ipjC 

= J2 (b A \H\c)\c)- J2 (b A \H\^c)\c), 

using the involutive nature of for the second sum. Now, using 2. from Lemma|9] 

(b A \H |^-c) = (-l) bA '^ c = (_!)& A -c+bi.- J+ i _ (_i)^- J+ i(b A |H|c) 

and thus 

■PlA^)= E (b A \H\c)\c)-(-l) b ^ £ (b A \H\c)\c). 



C<lf>jC 
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The conclusion in the case = 1 is now obvious. 

As for the case bL-j+i = 1, we see, using item 1. from Lemma|9j that 



« '- J+lb) >= X;<(^_ J+1 &) A |ff|c)|c> = -c| c > 



= ^(-l)^ (bA) ' c |c)= ^(-l) bA ' c +^| C ) 



£ (fo A |ff|c)|c) + £ (b A |if|c)|c) =-^1^). □ 



c:Cj— 



Corollary 10 /fw consider the special case where ctb = far all b £ M L , except aoo...o = cto and 
CK10...00 — a X> aR d where fii = ... = /3_l-i = 1 o«<^ Pl = P, then the determinant of the Theorem 
simplifies to the product III • II2 • II3 • II4 of the following four terms: 



Ii 1 = Yl (a + ax - 2k)( L ^) IT 2 = JJ (a + c*i - /3 - 2k + l)(^-i) 

0<2k<L 0<2fe-l<i 

n 3 = Yl (a - a 1 - 2k + n 4 = J} (a -ai -/?-2/c)( 1 ^ 1 ) 



0<2fc-l<i 0<2fc<i 

For the proof note that each b <E B L we get as the contribution from Al 
A b A • a. = a c ( fo A | H | c) 

= a (b A I i? I 00 . . . 00) + a x (b A \ H | 10 . . . 00) 

= a (-l) b - A - 00 - 00 + a,(-l) 6 - A U, - ° = a + a, 

because A • 00 ... 00 = 11 ... 11 and then b ■ 11 ... 11 = ||6|| mod 2, where ||6| denotes the Hamming 
weight of b and where we have used the fact that A is a symmetric matrix. Thus the 2 L eigenvalues are 

a + (-l) l|f, "a 1 -(3 zev -b (beM L ). 

Now there are four cases to consider: 

1. II&H is even and b\ = 0: this gives eigenvalues a$ + a\ — \\b 2 b 3 . . . fr^| and since bi does not 
contribute to ||6|| the vector 6263 . . . fri must have even weight 2k. There are ( i 2 ^, 1 ) possibilities 
which account for Hi. 

2. is even and b\ = 1: this gives eigenvalues «o + a\ — (3 — \\b2b3 . . . fri|| and since b\ does 
contribute to ||b|| the vector b 2 b 3 . . . b^ must have odd weight 2k — 1. There are (^tji) possibilities 
which account for II2 . 

3. is odd and b\ — 0: this gives eigenvalues ao — a± — ||i>2&3 ■ • - frill and since b\ does not 
contribute to ||b|| the vector 62&3 . . . fr_L must have odd weight 2k — 1. There are possibilities 
which account for II3. 
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4. is odd and b\ = 1: this gives eigenvalues ao — ct\ — (3 — \\b2b3 . . . bi\\ and since b\ does 
contribute to ||6|| the vector 6263 . . . 6l must have even weight 2k. There are ( L ^j}) possibilities 
which account for II4. 

Corollary 11 Setting now ao — —a and a\ — a, i.e., specializing as in Lemma^ gives for detM^ a 
product Tl[ ■ Tl' 2 ■ 1T3 • II4 of the following four terms 

n; = n (-2fc)( i - 1 ) 

0<2fc<L 

n 3 = J| {-2a -2k + 

0<2fc-KL 

which are precisely the 2 L distinct eigenvalues of the original Conjecture. 

5 Concluding remarks 

We have been able to solve Ayyer and Mallick's conjecture about the eigenvalues of the asymmetric 
annihilation process by embedding it into a more general model and using an orthogonal transform which 
makes the transition matrix upper triangular. In sharp contrast to the original problem, the general situation 
with parameters a b and bj (which may be given a "physical" interpretation using ( b | Al \ c) — a&e c 
and <[TJ,([2j>) is easier to handle because it is not degenerate: all "symbolic" eigenvalues are simple. Our 
proof does not seem to explain the maximum amount of degeneracy, as stated in Conjecture [T] 

On the other hand, we mention that the result of Corollary[3]about the partition function can be extended 
to the more general model. Again, in contrast to the inductive approach of Ayyer and Mallick in J2), as 
outlined in Sec. [3] we can solve this problem directly by transforming it orthogonally into the basis where 
it shows its triangular structure. 

We start by remarking that the columns sums of the extended model are constant a c , though not 
zero. This implies that ( 1^ | is the unique left eigenvector with eigenvalue a = a c of Mi(a, /3). The 
right eigenvector | x ) with the same eigenvalue corresponds to the steady state distribution of the original 
problem. Then | y ) = Hl \ x ) satisfies Ml (a, /3)| y ) = a\ y), where Ml = Hi ■ Ml(o., f3) ■ Hl is the 
matrix seen in the Wl -basis. This triangular system for y is written explicitly as 

(A* a +/3 rev -b) 2/5= (beB L ) 

j:b 3 =l 

where 

X* b = a-X b = 2 fl c (beM L ). 

c : b-c—1 

Note that the sum on the right only contains terms y c where c = ijjjb < b. For b = 00 ... the equation 
is void, so we may put J/00...0 = 1- Since the polynomials A^ A + /3 rev • b are mutually coprime, this shows 
by induction that the denominator of the rational normal form of y b is the product of all polynomials 
A* A + /3 rev ■ c, where c runs over the binary vectors that can be obtained from b by successive application 
of decreasing ^ -transformations. Consequently, the product of linear polynomials 



n 2 = [J (-/3-2fc + l)(^- 1 i) 

0<2fc-KL 

n 4 = Yl (-2ai -^-2/c)( L ^ 1 ) 



0<2fe<L 
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is the least common multiple of the denominators of the j/j,. This property is invariant under the Hadamard 
transform, so it applies also to the coefficients of | x ) = Hl \y). But 

(l L \x) = (l L \H L \y) = 2 L ' 2 (100 . . . 00|y) = 2 L ' 2 y o...o = 2 L/2 , 

so 2~ L / 2 | x) is already normalized and can be seen as the "symbolic" stationary distribution in the 
generalized model. What we have shown is: 

Theorem 12 Z(a, (3) is the partition function related to Mi, (a, 0). 

We conclude by remarking that the specialization as in Lemma [3] and Corollary [TT| brings us back to 
Corollary 3 This is not completely obvious, since the expression in Corollary 3 has only (^J 1 ) factors, 



whereas inTheorem 
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Jthere are 2 L — 1 factors. What happens is that upon specialization the requirements 
for least common multiples and greatest common divisors change. Taking this into account one finds that 
from the general expression for Z(a, f3) only the ( L ^ 1 ) terms where b <E M L with ||b| = 1 or 2 contribute 
- and this is precisely the statement of Corollary|3] 
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